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Abstract 



We present a simple technique for generating new solutions of Ein- 
stein's equations using such function transformations that leave the 
' field equations in the Ernst form. In this context we recover all the 

f^) . known covariant transformations of Ernst equations and we find the 

' role of the analytic ones. Finally we obtain a new asymptotically flat 

U" 1 solution starting from the Kerr solutions. 

U 

1 Introduction 

^ ■ An interesting class of solutions in general relativity are those stationary 

and axisymmetric. In particular the ones which are asymptotically flat are 
candidate for the description of the exterior gravitational field of a uni- 
formly rotating star. The only known stationary axisymmetric asymptot- 
ically flat solutions are the Kerr solutions pQ, the generalized Tomimatsu 
solutions [2j and the Papapetrou solutions One would like, therefore, 
to obtain more such solutions. Ehlers first [I] showed how it is possible to 
construct new stationary exterior solutions and stationary interior solutions 
of the Einstein's field equation starting from static exterior solutions by 
applying certain conformal transformations to auxiliary metrics defined on 
three-dimensional manifolds in space-time. These are one-parameter fam- 
ily of solutions. Geroch [5] has shown that one can, in fact, obtain an 
infinite-parameter family. Further Xanthopoulos [H] has given a technique 
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for generating from any one-parameter family of vacuum solutions of Ein- 
stein's equation a two-parameter family. In this paper we introduce a simple 
but very powerful method for generating an infinite tower of solutions. We 
start from a known solution using function transformations on the complex 
potential e = / + which leaves the Ernst equations covariants in form. 
The plane of this paper is the following. In section 2 we make an overview 
on the Ehlers method for vacuum solutions in relation with the Papapetrou 
gauge for stationary axisymmetric space-time and we find a functional re- 
lation between / and <£. In sections 3-5 we derive the basic equations and 
in section 7 these are appropriately reduced. In sections 6,8-9 we integrate 
the basic equations in simple cases. Finally, in section 10 we derive a new 
asymptotically flat solution starting from the Kerr solutions. 



2 Ehlers method for the vacuum 

The Ehlers method [3] makes possible to construct new stationary exterior 
and interior solutions of the field equations of general relativity starting from 
static solutions. Ehlers method for the vacuum affirms that if u = Uidx 1 is 
a linear differential form and we take a solution of the static form 

G = e~ 2U H - e 2U dt 2 , V 2 U = 0, (1) 

where H = hikdx % dx k is a quadratic differential form, is then possible to 
obtain solutions of the form 

G = cosh(2C/)# - cosh" 1 (2[/) (eft - u) 2 , (2) 

provided that the following condition is satisfied 

VjkiU' 1 = Ufa k ], (3) 

where denotes the usual totally skew-symmetric tensor with components 
0, ±y/h , h = \hik\ and "," denotes the ordinary derivate . In relation to 
equation © for the vacuum axisymmetric case we have, in Papapetrou 
gauge [H], tii = U2 = , U3 = u, with 

H = e^Kdx 1 ) 2 + (dx 2 ) 2 } + p 2 {dx 3 f. (4) 

Starting from the Ernst equations [Zj with / _1 = cosh(2?7) and with the 
equation for lo given by 

UAA ~ u A h -j^a/a = , A = 1, 2, (5) 

p f 
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it is easy to show that the solution of © and © is given by 



ux = 2pU 2 , lo 2 = -2 P Ul (6) 

From equation © it follows that V 2 U = 0. For the metric function / we 
have 

2U A = j± ■ (7) 

From the Ernst potential $ given by 



it follows that 



and therefore 



$ A = P^— (9) 



$ = V 7 ! - / 2 + c (10) 
The constant c is inessential. Finally we find 

/ 2 + $ 2 = l. (11) 

The condition 1)11(1 caracterizes the vacuum stationary axisymmetric solu- 
tions achievable with the Ehlers method compatible with the Papapetrou 
gauge. The condition (fTT|) also caracterizes the Papapetrou class of solutions 
[3]. In Appendix A we give some solutions derived from the relation (|llj) . In 
section 8 we see how a "gauge" version of the condition (|11|) emerges in the 
context of our method, that can be therefore considered a generalization of 
the Ehlers method for the vacuum case. It is easy to show that the condition 
(|11|) leads to solutions which are the imaginary part of the Ernst £ function 
with e = f + %<$> = 

3 Some basic facts 

If we take the Ernst equations [7j 

/V 2 / + ^ A -fl = 0, 

/V 2 $ - 2f A $ A = (12) 
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and if we perform, based on (|12|) . a transformation in the new variables a 
and b such that 

aV 2 a + b\ - a\ = 0, 

aV 2 b - 2a A b A = 0, (13) 

then the transformation / = f(a,b), = <3?(a, b) maps old solutions / 
and <& in the new solutions a and b, where a and b are the new / and 
<]?, respectively. Viceversa, if a and 6 are old solutions, / and $ are new 
solutions. These transformations provide significantly new solutions derived 
from the circumstance that an integration is involved in passing from <3? to 
the metric functions. 

4 Formulation of the problem 

Beginning from (|12[) we pose / = /(a, b) and $ = $(a, 6), from which, after 
a simple but very tedious algebra, (where index A = 1,2 corresponds to the 
spatial derivates with respect to x 1 and x 2 and index a and 6 denote the 
derivate with respect to the variables a and b) it follows 

/ JV 2 a + L<4 + Mb\ + 2Na A b A = 0, (14) 
/ JV 2 6 + Lai + Mb 2 A + 2J\TaA&j4 = 0, (15) 

where J = / a $ b - f b ® a and 

£ = fifaa^b ~ *aafb) + " /a) + 2/aA$a, (16) 

M = /(/ 66 $ 6 - $ M / b ) + + fb), 

N = f(f ab $ b - <5> ab f b ) + $ a (/ fc 2 + $ 2 ), 

L = -fUaa^a ~ <Wa) ~ *a(/a + *«), 

M = -/(/»*„ " ®bbfa) ~ <M<^ ~ fi) ~ 2f a f b * b , 

N = -f(fab$a ~ <Wa) ~ $b{fl + 

Note that 

fJa = (fa®b - fb$a) a = f(faa$b ~ fb^aa) ~ f(fab®a ~ ®abfa), 
but from equation (|16|) we find 

l+ ^ /3 (?1- <i7) 
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In the same way from fj b = f(f a ^b — fb^a) b we obtain 

J 



N + M = f a [ ■ (18) 



5 Basic equations 

We define 



X = j f a , Y = jf b , £ = j$ a , V = j$ b (19) 



and then 



L = f[r } X a -Yt a + r } e + XYZ], 
M = f[7 1 Y b -Y Vb + r ] ( V 2 + Y 2 )}, 

N = f —[q{X b + Y a ) - F(6 + Va) + Z(V 2 + Y 2 ) + r)(XY + tf)], 

L = -f[iX a -Xi a + i{x 2 + e)l 
M = -f m - Xrj b + 7 1 (XY + £77)] , 

N = -^[£(X b + Y a ) - X{£ h + ria) + 7!(X 2 + f) + i{XY + £77)]. 

From equation (|12|1 taking / = f(a,b) and $ = <&(a, 6) we obtain equations 
(HU-CSJ with 

J = f(X V - Y£) ? 0. (20) 

For the covariance to arise and consequently for equations p4 |) -([15 |) to take 
the form (|13j). first we must have 

L = M = N = 0. (21) 

Besides we have equations (|17 j) -([18 j) with the condition 

-L = M = -N (22) 

and therefore equation (|17|) becomes 

2N '=f i (P) ■ (23) 



with 

2 v/ 



L = -M = V( Tf ) . (24) 
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While for equation (JT8|) we get 



N = r 



(25) 



and then 



/ 3 72 v '« + — ( 72 (4-^a)+/ 3 ( 71 ) «a6a = 0, (26) 



J 



P + 2 [pj a ^ 1 J V/ 2 / b 



±X72 h _flfl 

p 2 V/V6^ ~^ 1 J v/v; 



/ 3 ^6-V(^) (4-&i)+/ 3 (4) <M^ = (27) 



and finally 



J 



J 



1 



(28) 



The equations (|26|) - (|27|) are the same as equations (|13|) and therefore covari- 
ance follows. 



6 Case with £ = aX 

At the first stage we consider the equation L = 

£X a -X£ a + £(X 2 + £ 2 ) = 0. 

The equation (|29j) permits to take 

i = aX 

and from equation (|29|) we have, after an integration 



Q: 



e = ±1. 



(29) 



(30) 



(31) 



Where B = B(b) is an integration constant. In the same way, taking equa- 
tion (|3"Tj) into we obtain 



5 - $ = e\/S 2 - / 2 , 



(32) 



where i? also depends only by 6. Deriving equation (|32|) with respect to b 
we obtain 



V 



1 

7 



-B& — eBB b 



1 



(33) 
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The four variable X, Y, £, rj are then related by 

f n n f 

Z = aX = -^X,r ] = - + aY = - + 1 -L-Y, (34) 

where 

n = B b - BB b -^ — , n' = BB b - eBB b - B b $. (35) 
B - <f> 

Beyond the six equations L = -M = N = / 3 /2(J// 2 ) a and L = M = N = 
we have condition (|2U)) 

Xr] - Y£ = A = -i (36) 

but from equations (j34j) we have 

A U Af 

from which it follows 
1 

~7 



B b --B b ^W^P = (38) 

B J a 



We set A* = l/a and /i(6) = -B*. In this case we can obtain / as a function 
of A* , B*, B b and B. With a little algebra from (jSHJ) we find 

fif 2 -2uf + 9 = 0, (39) 

where 

= (A* + B*f + j^Bl , i/ = B 6 (A* + B*) , 9 = Bl - B 2 b . (40) 
Finally from the equation ()39|) we obtain 

f = u±V^EE. (4i) 

Note that we have obtained the equation (|4*Tj) using only the equation L = 0. 
We must therefore check this solutions with the others equations. In the next 
section we reduce the relevant equations. 
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7 Reduction of the relevants equations 

We start from the six equations quoted above 

n(X a + e) - Y{i a - XO = ~A a , (42) 

r,(Y b + r, 2 )-Y( Vb -Y V ) = ~A a , (43) 
V [(X b + Y a ) + 2^] - Y[(£ b + n a ) - (Xr, + Y£)] = 0, (44) 

S(x a + e)-Y(£ a -xt;) = o, (45) 

£(Y b + if ) - X( Vb - Yn) = 0, (46) 
£[{X b + Y a ) + 2^] - X[i b + Va - (Xr, + Y$] = -A*, (47) 

with condition (|36|) . But f(a,b), $(a,6) are not arbitrary. The following 
condition must be imposed 

fba = fab , <S>ab = $b a , (48) 

from which it follows 

X b = Y a , & = r) a + A. (49) 

Besides, differentiating condition (|3fij) we have 

r/X a - y£a - [Cy a - Xr/J = A a , (50) 
r/X 6 - Y& = en - Xt7 6 . (51) 

Using equations (|49|). equation (|H|) becomes 

2[(7?y 6 - + e?? 2 + = o (52) 

and thanks to equation (|51|). equation (|o^|) is the same as (|46j). Similarly 
equation (fl9|) placed into (|4T|) leads to 

2[£Y a - X Va ] + 2f V -XA + X(X V + YO = 
= 2[tY a - X Va + 2f V + 2XY£] = 
= 2[r/X a - Yi a -A a + V e + XYi] 

r/(x a + e 2 )-y(ea-xo = ^ a . (53) 
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But the equation is the same as the (|4*2*|) . Therefore in place of equations 
1(12 J) -P7 | we have 



V (X a + f) - Y{i a - XO = l -A a , (54) 

r 1 (Y b + r) 2 )-Y( Vb -Yr l ) = --A a , (55) 

i{X a + f ) ~ X(U - XQ = 0, (56) 
£(Y fe + rf) - X( Vb - Yr,) = 0, (57) 

with the condition (|36|), Thus the equations for the variable a and b are 
the Ernst ones if and only if equations (|54 |l -l|57 |l are satisfied. Note that 
this system is a first order system in the derivate a and 6, while the Ernst 
equations are a second order system in the coordinate derivates. In the next 
section we conclude the integration in the case with £ = aX. 



8 Final solution in the case with £ = aX 



At this stage we have used equation (|56f) only. From the second equation of 
(jE3J, with the help of (JHOJ),®, we get 



Y = 7] 



A 

T 



(58) 



from which, after a derivation, we obtain 

„2\ 



i{Y b + r,') - X( Vb -Y V )=Z 



f 



A 



If equation 1)57(1 is true, then it follows (£ ^ 0) 



/ 



A' 

o 6 



(59) 



(60) 



We suppose first that equation ((6fl|) is verified, then equations (|54 |) -(|55 )1 

become 



X 



1 



Vj ~ Y ) (6. - = T^a 

- y ) (.Vb - Yr)) = ~A a . 



(61) 
(62) 
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Equation (|6*Tj) gives 



B = B (b). 



(63) 



In relation to equation (|41jl and thanks to equation (|35|) . (|37|) after a little 
algebra, we find 



IT = -efJB 2 (A* + B*Y 



B : 



from which we get 



with the only solution given by 

B* = . 



, B = ± 



1 

B' 



(64) 



(65) 



(66) 



We must consider the last two equations. At this stage we have the equations 



f + $ 2 = 2B<S>e , / 



B b 



1 + 3 ' 



(67) 



Note that the condition (|6*Tj) for / and can be considered a "gauge" version 
of the Ehlers condition (jllj) . The remaining equations are the equation (|6()|) 
and (|62|). For equation (|6*2*)) we have 



i *7 

ln 7 



2 A 



From the relation ()63|) . after an integration, we find 

77 + g(q) 



(68) 



(69) 



where G = J ^db and H = H(a) is an arbitrary integration function. Plac- 
ing equation (|69|) into (|6Uj) we obtain 



[G + H]B b 



also, from (j63j) . we have 



^ [/ 6 £ - /5 6 ] . 



(70) 



(71) 
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With the expression for / given by (|67|). equation ffTTTf) gives 
[G + H}B b = 



+ 



a 2 
2^ 2 



1 r R2 r r \ B 6 , B bb B l 

-{B b -BB bb )- W2 +^- 



(72) 



Equation (|72() has two solutions. The first is the one with H = 0, from 
which it follows that fi is a constant, i.e. / = 0, and therefore it is not 
acceptable. The other solution is given by 

GB b = ^(l-^) ,H = ea\B b = 1 -^ T ^l, (73) 

from which it follows 

-l = 2e6 + c. (74) 

The constant c is inessential. Therefore we can choose c = 0; in this case 
G = -eb 2 . We get 

f = -T^rjt > 72 *b = -A + e ' > e ' = ±L ( 75 ) 
From the first equation of (|fi7|) we have 

$ 2 + ef + / 2 = 0^ 
o 

»■ a 2 



$ 1 = 5 i $2 = -77^3 PTT- (76) 

a 2 + b 2 b{a 2 + b 2 ) ' 



Taking now the solution for $i and placing into (|75|) we have 

b 2 - a 2 _ -b 2 e + e'a 2 

9 9 ' 



(77) 



from which taking e = — 1 and e' = —1, it follows that the equation l)75j) 
is verified. It is easy to show that the equation £ = (l//)$ a is verified. In 
the same way we can see that the solution $2 does not verify the condition 
1)75(1 . Then we have the final solution 

/ = 9 a To * = o-^T?- (78) 

a 2 + b 2 a 2 + b 2 

The solution ((78|) represents the inverse transformation, i.e. it says that if 
e = / + «<!> is a solution, also 1/e is a solution. Although this is a known 
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transformation, it has been obtained with a new method. Our starting 
point is equation (|56|) with the assumption £ = aX. If we consider a more 
general situation, i.e. £ = F(X,Y,r/), we can obtain more general solu- 
tions. Therefore a generalization of condition (|67|) appears. In this sense 
our method can be considered a generalization of the Ehlers method. In 
the case with £ = it is easy to show that the only solution is the trans- 
formation f = a , = b + const, which is the identity transformation. It 
is reassuring that if we consider the most simple ansatz we obtain the most 
simple solution. The search of new transformations is an open question, but 
our method gives a tool in such direction. In the next section we consider 
the analytic case, i.e. F = —Y. 

9 Analytic case 

In the analytic case we have for the variable X,rj,Y,£ the conditions 

X = n , y = -e, (79) 



with 



a 

The main equations in this case are 



X V-Y£= 72 ^X 2 + Y 2 = -i (80) 



X(X a + Y z )-Y(XY -Y a ) = - — , (81) 

cr 

X(Y b + X 2 )-Y(X b -XY) = \, (82) 

a 6 

XY a - YX a - Y(X 2 + Y 2 ) = 0, (83) 

-Y(Y b + X 2 ) = X(X b - XY). (84) 

At first we check the consistency of our equations. From equation (|84[) we 
have 

X b -XY = -^(Y b + X 2 ), (85) 
which taken into (|82[) and using equation ()8(J|) gives 

^2 X 



Y b + X z = —. (86) 



a 



Equation (|81|) is expression (|8()[) differentiated with respect to a. From 
equation (fSUl) we have 

X = -Vl~a 2 Y 2 . (87) 
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In the expression (|87|1 we have choosen the plus sign. We have the conditions 
Q49JI from which it follows 



n 



From equation (|83j) we obtain 



X„ 



X, 
Y 



1 



(89) 



and then expression Q88JI becomes 

n = - Jy a - (90) 
Differentiating expression (|87j) with respect to a, thanks to (|9"U|) . we have 



^a = -4r\/l-a2y2 



y 2 



a 2 Y 2 Y h 



Vi - a 2 r 2 i-a 2 Y 2 ' 



+ 



(91) 



which taken into (|88|) , after a little algebra, it gives equation (|86[) . Therefore 
our equations are consistent. From equation (|86|). after an integration, we 
get 

Y 2 [a 2 + (b + H(a)) 2 ] + 2Y[b + ff(o)] = 0, (92) 

where H = H(a) is an integration constant. A simple solution is Y = 
0, which leads to the identity transformation, a obviously trivial analytic 
transformation. The other solution is 



Y 



2[b + H(a)} 
[a 2 + (b + H{a)f] 

and then from f b /f = Y we obtain 

p G(a) 

f 



(93) 



[a 2 + (6 + #) 2 ] ' 

We must finally consider the condition (|8U|) with X = f a /f. We obtain 



(94) 



G 



+ 



4aG a 



[a 2 + (6 + tf) 2 ] [a 2 + (6 + if) 2 ] 



4ff Q (6 + ff) 
[a 2 + (6 + ff) 5 



g a (6 + H) 
[a 2 + (b + H) 2 } 



+ G a + a 



(95) 
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Because H and G depend on the variable a, we conclude that the only 
solution of the equation (j9"5|) is obtained posing G = In a and H = c, where 
c is a constant. If we choose c = we obtain again the transformation 

/ = 7fT» < 96) 

and for the analyticity condition we must have 

(97) 

which represents again the inverse transformation. It is a simple matter to 
verify that choosing the minus sign in expression Q87|) we can never obtain 
an analytic solution. We set now in equation l)94|) H = — c and therefore we 
obtain 

f = a 2 + b 2 +c 2_2bc ^ 

With simple manipulations (if a and b are solutions, also ca and cb are 
solutions) we see that the solution (|98|) can be written in the following way 



which for $ it gives 



f c 2 a 2 + c 2 h 2 + l _ 2d) ' (") 



_ b-cb 2 -ca 2 

l + c 2 a 2 + c 2 b 2 -2cb' ( ' 



The solutions (|9*9*)) - H1UUI) say that if e = a + «6 is a solution, also it is a 
solution the transformation 

(101) 



1 + ice 



The solution (|1()1[) represents the so called Ehlers transformation. We can fi- 
nally conclude with the claim that the only analytic, not trivial, transforma- 
tions which allow the covariance of the Ernst equations, are the inverse trans- 
formations and the Ehlers transformations. Note that we can not obtain an 
infinite tower of solutions starting from the inverse transformation, because 
the inverse transformation of the inverse one is the identity transforma- 
tion. For the Ehlers transformation we have that, after n-transformations, 
e — > tt^ — • We can start from two arbitrary transformations to obtain 
infinite combinations that allow infinite covariant transformations. In the 
next section we apply the inverse transformation to the Kerr solutions 
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10 Function transformations for the Kerr metric 



We consider for simplicity the case of the extreme Kerr solution JT] expressed 
in spherical coordinates with r = e v . We start from the line element 



ds 2 = f-^e^idv 2 + cM 2 ) + p 2 dip 2 } - f(dt - tudip) 2 . 
We have for this solution 

p"e" 



O — fKerr 
b = <&Kerr 



(pe v + l) 2 + cos 2 i?' 

2 cos $ 
(pe v + l) 2 + cos 2 tf' 



(102) 

(103) 
(104) 



Performing the inverse transformation we obtain for the new / and the 
expression 



[p 2 e 2v - sin 2 tf] {pe v + l) 2 + cos 2 ?? 



(p 2 e 2v — sin 2 i?) 2 + 4cos 2 t? 



2 cos $ 



{pe v + l) 2 + cos 2 i? 



(p 2 e 2v — sin 2 i?) 2 + 4cos 2 $ 



(105) 



We considere now the asymptotical behaviour of the solution (|1U5|) . First 
we note that for v — * oo / — * 1 and —* 0. Besides we have 



2sini? . _ s 4 cos i? 



p 2 e 2v 



from which, thanks to (JBJ), we get 

2 



a; 



p 2 e 2v 



sm z i}e- v + 0(e~ 2v ) 



(106) 



(107) 



and therefore for v — > oo uj — > 0. From the Ernst equations [Jj it follows 
that 

^sin^e-^ + CKe- 3 ^). (108) 



7 ~ v 



2p 2 



Setting r = e v we obtain 
2 

pr 



2 



f ~ 1 H , cj ~ — jSm 2 $ ) 7 — mr 



2p 2 r 2 



(109) 
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Therefore the solution ()1(J5|) . obtained performing the inverse transforma- 
tion to the extreme Kerr solution, is asymptotically flat. Besides, thanks 
to expansion (|1U9|) . the solution Q1U5[) represents a metric with a spinning 
source if and only if p = —1/m, where m is the mass of the source. Similar 
arguments follow for the other Kerr solutions. The application of the Ehlers 
transformations to the Kerr metric does not give asymptotically flat solu- 
tions. We conclude with a remark on the Ehlers transformations. Consider 
now the transformation <|1U1|) . We consider the case with c = 1. Performing 
this transformation to the class of solution with f 2 + <I> 2 = 1, we obtain 

f ' (110) 



2(1-$) 2 

Because the constant c is an inessential constant (for c ^= 1, if e is a solution 
also | is a solution) we can say that performing a Ehlers transformation to 
the class of solutions given by the condition (jlljl (c = 1 or e — * ~) we obtain 
a static solution (up to a constant factor which can be set to 0). 



11 Conclusions 

In this paper we introduce a new, simple but very powerful method for 
generating a tower of solutions starting from known ones using the covariance 
of the Ernst equations for generic metric functions transformations. The 
basic equations are derived and reduced. We have considered particular 
cases obtaining all the transformations known in literature. Finally, we 
obtain an asymptotically flat solution starting from the Kerr one. 
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Appendix A 

In relation to we consider first the cylindrical solution $ = &(p) and 
therefore u = lo(z). This ansatz gives / = ^f\ P with p real constant and 
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integrating the field equations we have 

ds 2 = ^-(pP + p -P)(dp 2 + dz 2 ) + 

+ y (p* + P' V W ~ -^f^{dt-pzdyf. (Ill) 

This is a cylindrical solution which is an anti-Machian solution. 
It is possible to give a non cylindrical solution starting from the equation 
fj 1 1 1 1) using parabolic coordinates p = 2Xp, z = X 2 — p? . The Laplacian in 
these coordinates is 



V 2 = ^ + AA + ^ + ^, (112) 

A p 



from which it follows that if the metric functions are indipendent from A 
(or p) the Laplacian (|112j) is isometric to the cylindrical one. Therefore if 
/ = /(A) we obtain, similarly to the above cylindrical solution, 

ds 2 = ^±^±^l(dX 2 + dp 2 ) + 

2xp(^y 4 

2p 2 (l + X 2 P) , 2 2XP 2 2 

+ T^Z2 dl P ~ I T X 2 P ( dt ~ d( P) ■ (113) 

Equation (|113|) has a complicated behaviour when expressed in cylindrical 
coordinates. If we change A with p in the line element ()113|) we have also 
a not cylindrical solution. By using spherical coordinates U = —k/r is a 
solution of V 2 U = 0. Therefore we have the spherical solution 

ds 2 = e " iL T2 L - i cosh(2A;/r)[dr 2 + r 2 d6 2 } + r 2 sm 2 6 cosh(2k/r)dip 2 - 

(dt - 2k cos 6 "dip) 2 . (114) 



cosh(2/c/r 

Finally, taking U = pz we obtain the solution 

ds 2 = cosh(2pz)e~ p2p2 [dp 2 + dz 2 } + p 2 cosh(2pz)d(p 2 - 

^r—(dt-pp 2 dpf. (115) 

cosh(2pz) 

Note that taking z = in (|115[) we recover the Van Stokum solutions jSj in 
a co-moving reference frame. 
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